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Abstract 



O 

CN We employ the Yukawa coupling unification condition, y t = y b = y T at M GU T5 

inspired by supersymmetric SO(10) models, to estimate the lightest Higgs boson 

^ mass as well as masses of the associated squarks and gluino. We employ non- universal 

^vq soft masses, dictated by SO(10) symmetry, for the gauginos. Furthermore, the soft 

masses for the two scalar Higgs doublets are set equal at Mqut, and in some examples 

, , these are equal to the soft masses for scalars in the matter multiplets. For \i > 0, 

M 2 > 0, where M 2 denotes the SU (2) gaugino mass, essentially perfect t-b-r Yukawa 
unification is possible, and it predicts a Higgs mass of 122 - 124 GeV with a theoretical 
uncertainty of about ±3 GeV. The corresponding gluino and the first two family 

,-£< squarks have masses > 3 TeV. We present some LHC testable benchmark points 

which also show the presence of neutralino-stau coannihilation in this scenario. The 
well-known MSSM parameter tan (3 ~ 47. 
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1 Introduction 



Supersymmetric (SUSY) 5*0(10) grand unified theory (GUT), in contrast to its non- 
SUSY version, yields third family (t-b-r) Yukawa unification via the unique renor- 
malizable Yukawa coupling 16 • 16 ■ 10, if the Higgs 10-plet is assumed to contain 
the two Higgs doublets H u and of the minimal supersymmetric standard model 
(MSSM) [1]. The matter 16-plet contains the 15 chiral superfields of MSSM as well 
as the right handed neutrino superfield. The implications of this Yukawa unification 
condition at Mqut ~ 2 x 10 16 GeV have been extensively explored over the years 
[1, 2]. In 50(10) Yukawa unification with > and universal gaugino masses, the 
gluino is the lightest colored sparticle [3, 4], which is now being tested [5] at the Large 
Hadron Collider (LHC). The squarks and sleptons, especially those from the first two 
families, turn out to have masses in the multi-TeV range. Moreover, it is argued in 
[3, 4], based on the results of publicly available codes like Isajet [16] and Softsusy 
[6], that the lightest neutralino is not a viable cold dark matter candidate in 5*0(10) 
Yukawa unification with fi > and universal gaugino masses at Mqut- On the other 
hand, in Ref.[3] it is shown that neutralino DM is allowed through h-resonance for 
Yukawa unification at the 10%level. 

Spurred by these developments we have investigated t-b-r Yukawa unification [4, 
7, 8] in the framework of SUSY SU(A) e x SU(2) L x SU(2) R [9] (4-2-2, for short). The 
4-2-2 structure allows us to consider non-universal gaugino masses while preserving 
Yukawa unification. An important conclusion reached in [4, 7] is that with same 
sign non-universal gaugino soft terms, Yukawa unification in 4-2-2 is compatible with 
neutralino dark matter and gluino co-annihilation [4, 5, 7, 10] being a unique dark 
matter scenario for /i > 0. 

By considering opposite sign gauginos with fi < 0, M 2 < 0, M 3 > where fi is the 
coefficient of the bilinear Higgs mixing term, M 2 and M 3 are the soft supersymmetry 
breaking (SSB) gaugino mass terms corresponding to SU{2) l and SU(3) C respectively, 
it is shown in [8] that Yukawa coupling unification consistent with the experimental 
constraints can be implemented in 4-2-2. With /i < and opposite sign gauginos, 
Yukawa coupling unification is achieved for m 16 > 300 GeV, as opposed to m 16 > 
8 TeV for the case of same sign gauginos. The finite corrections to the b-quark mass 
play an important role here [8]. By considering gauginos with M 2 < 0, M3 > and 
fi < 0, we can obtain the correct sign for the desired contribution to (g — 2) M [11]. 
This enables us to simultaneously satisfy the requirements of t-b-r Yukawa unification 
in 4-2-2, neutralino dark matter and (g — 2)^, as well as a variety of other bounds. 

Encouraged by the abundance of solutions and co-annihilation channels available 
in the case of Yukawa unified 4-2-2 with M 2 < and fi < 0, it seems natural to 
explore Yukawa unification in SO(10) GUT (with M 2 < and fi < 0). It has been 
pointed out [12] that non-universal MSSM gaugino masses at Mgut can arise from 
non-singlet F-terms, compatible with the underlying GUT symmetry such as SU(5) 
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and SO(10). The SSB gaugino masses in supergravity [13] can arise, say, from the 
following dimension five operator: 



2M P 



A a A 6 + c.c. 



(1) 



Here A a is the two-component gaugino field, F denotes the F-component of the 
field which breaks SUSY, the indices a, b run over the adjoint representation of the 
gauge group, and Mp = 2.4 x 10 18 GeV is the reduced Planck mass. The resulting 
gaugino mass matrix is (F ab ) /Mp where the supersymmetry breaking parameter (F ab ) 
transforms as a singlet under the MSSM gauge group SU(3) C x SU (2) L xU (l)y. The 
F ab fields belong to an irreducible representation in the symmetric part of the direct 
product of the adjoint representation of the unified group. This is a supersymmetric 
generalization of operators considered a long time ago [14]. 
In SO (10), for example, 



If F transforms as a 54 or 210 dimensional representation of SO(10) [12], one obtains 
the following relation among the MSSM gaugino masses at Mgut : 



where Mi, M 2 , M 3 denote the gaugino masses of U(l), SU{2) L and SU(3) C respec- 
tively. The low energy implications of this relation have recently been investigated 
in [15] without imposing Yukawa unification. We consider the case with fi > and 
non-universal gaugino masses defined in Eq.(3) in this paper. In order to obtain the 
correct sign for the desired contribution to (g— 2) M , we set signs of the gaugino masses 
as Mi > 0, M 2 > and M 3 < 0. Somewhat to our surprise, we find that this class of 
t-b-r Yukawa unification models make a rather sharp prediction for the mass of the 
lightest SM-like Higgs boson. In addition, lower mass bounds on the masses of the 
squarks and gluino are obtained. 

Notice that in general, the soft terms such as the trilinear couplings and scalar 
(mass) 2 terms are not necessarily universal at Mgut- However, we can assume, 
consistent with SO(10) gauge symmetry, that the coefficients associated with terms 
that violate the SO(10)-invariant form are suitably small, except for the gaugino 
term in Eq.(l). We also assume that D-term contributions to the SSB term are much 
smaller compared with contributions from fields with non-zero auxiliary F-terms. 

The outline for the rest of the paper is as follows. In Section 2 we summarize 
the scanning procedure and the experimental constraints that we have employed. 
In Section 3 we discuss the important role of SUSY threshold corrections in t-b-r 
Yukawa unification. In Section 4 we discuss how radiative electroweak symmetry 



(45 x 45) 5 = 1 + 54 + 210 + 770 



(2) 



M 3 : M 2 : M x = 2 : 



3 : 



1, 



(3) 
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breaking (REWSB) is compatible in our model despite setting m Hu = m Hd at M GU t- 
In Section 5 we present our results, focusing in particular on the mass of the lightest 
(SM-like) CP even Higgs boson, as well as the masses of squarks and gluino. The 
table in this section lists some benchmark points which will be tested at the LHC. 
Our conclusions are summarized in Section 6. 

2 Phenomenological Constraints and Scanning Pro- 
cedure 

We employ the ISAJET 7.80 package [16] to perform random scans over the funda- 
mental parameter space. In this package, the weak scale values of gauge and third 
generation Yukawa couplings are evolved to Mqut via the MSSM renormalization 
group equations (RGEs) in the DR regularization scheme. We do not strictly enforce 
the unification condition g 3 = gi — g 2 at Mgut, since a few percent deviation from 
unification can be assigned to unknown GUT-scale threshold corrections [17]. The 
deviation between gi = g 2 and g% at Mgut is no worse than 3 — 4%. For simplicity we 
do not include the Dirac neutrino Yukawa coupling in the RGEs, whose contribution 
is expected to be small. 

The various boundary conditions are imposed at Mqut and all the SSB parame- 
ters, along with the gauge and Yukawa couplings, are evolved back to the weak scale 
Mz- In the evaluation of Yukawa couplings the SUSY threshold corrections [18] are 
taken into account at the common scale Msusy = ^ m t L m t R where m^ L and m^ R are 
the third generation left and right handed stop quarks. The entire parameter set 
is iteratively run between Mz and Mgut using the full 2-loop RGEs until a stable 
solution is obtained. To better account for leading-log corrections, one-loop step-beta 
functions are adopted for gauge and Yukawa couplings, and the SSB parameters m ; 
are extracted from RGEs at multiple scales m 8 = mi(mi). The RGE-improved 1-loop 
effective potential is minimized at Msusy, which effectively accounts for the leading 
2-loop corrections. Full 1-loop radiative corrections are incorporated for all sparticle 
masses. 

The approximate error of ±3 GeV in the ISAJET estimate of the Higgs mass 
largely arise from theoretical uncertainties in the calculation of the minimum of the 
scalar potential, and to a lesser extend from experimental uncertainties in the values 
for m t and a s . 

An important constraint comes from limits on the cosmological abundance of sta- 
ble charged particles [19]. This excludes regions in the parameter space where charged 
SUSY particles become the lightest supersymmetric particle (LSP). We accept only 
those solutions for which one of the neutralinos is the LSP and saturates the WMAP 
bound on relic dark matter abundance. 
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The MSSM Higgs doublets reside in the 10 dimensional representation of SO (10) 
and fermions of the third family belong to the 16 dimensional representation of 
SO(10), which implies Yukawa coupling unification at Mgut- 

We have performed random scans for the following parameter range: 

< mie < 5 TeV 
< m w < 5 TeV 
< M 1/2 < 2 TeV 
35 < tan (3 < 55 
-3 < A /m 16 < 3 (4) 

(5) 

Here m 16 is the universal SSB mass for MSSM sfermions, mio is the universal SSB 
mass term for up and down MSSM Higgs masses, M1/2 is the gaugino mass parameter, 
tan f3 is the ratio of the vacuum expectation values (VEVs) of the two MSSM Higgs 
doublets, Aq is the universal SSB trilinear scalar interaction (with corresponding 
Yukawa coupling factored out). We use m t = 173.1 GeV [20]. Note that our results 
are not too sensitive to one or two sigma variation in the value of m t [21]. We use 
iTib(mz) = 2.83 GeV which is hard-coded into ISAJET. 

In scanning the parameter space, we employ the Metropolis-Hastings algorithm 
as described in [22]. The data points collected all satisfy the requirement of REWSB, 
with the neutralino in each case being the LSP. After collecting the data, we impose 
the mass bounds on all the particles [19] and use the IsaTools package [23] to imple- 
ment the various phenomenological constraints. We successively apply the following 
experimental constraints on the data that we acquire from ISAJET: 



m h (lightest Higgs mass) > 114.4 GeV [24] 

BR(B S fi+fj,-) < 1.2 x 10~ 8 [25] 

2.85 x 10- 4 < BR(b -»■ si) < 4.24 x 10~ 4 (2a) [26] 

0.15 < B S^ t) ^ ssm < 2.41 (3a) 26 

— BR(B u ^tv t )sm — \ / L J 

Q CDM h 2 = 0.1123 ±0.0035 (5a) [27] 

0<A(#-2) M /2 <55.6xl0- 10 [11] 



Employing the boundary condition from Eq.(3) one can define the MSSM gaugino 
masses at Mqut i n terms of the mass parameter M1/2 : 

Mi = Mi /2 
M 2 = 3M 1/2 

M 3 = -2M 1/2 (6) 

Note that M 2 and M3 have opposite signs which, as we show, is important imple- 
menting Yukawa coupling unification to a high accuracy. 
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3 Threshold corrections and Yukawa unification 



The SUSY threshold corrections to the top, bottom and tau Yukawa couplings play 
a crucial role in t-b-r Yukawa coupling unification. In general, the bottom Yukawa 
coupling yf, can receive large threshold corrections, while the threshold corrections to 
y t are typically smaller [18]. The scale at which Yukawa coupling unification occurs 
is identified with Mqut, the scale of gauge coupling unification. Consider first the 
case ^(Mgut) ~ 2/t(^gut)- The SUSY correction to the tau lepton mass is given by 
Sm T = v cos f35y T . For the large tan values of interest here, there is sufficient freedom 
in the choice of Sy T to achieve y t f» y T at Mqut- This freedom stems from the fact 
that cos/3 ~ 1/ tan/3 for large tan/3, and so we may choose an appropriate Sy T and 
tan/3 to give us both the correct r lepton mass and y t ~ y T . The SUSY contribution 
to Syb has to be carefully monitored in order to achieve Yukawa coupling unification 
y t (M GVT ) « y b (M GVT ) « y T (M GVT ). 

We choose the sign of Syi (i = t, b, r) from the perspective of evolving yi from 
M GU t to Mi- With this choice, Syb must receive a negative contribution (—0.27 < 
fyb/yb i$ —0.15) in order to realize Yukawa coupling unification. This is a narrow 
interval considering the full range of —0.4 < Syb/yb < 0.4. The dominant contribution 
to Syb arises from the finite gluino and chargino loop corrections, and in our sign 
convention, it is approximately given by [18] 



r finite ^ 9l /imgtan/3 y\ /M^tan/3 
yb ~ 12vr 2 ml 32tt 2 m 2 ' {) 



where g% is the strong gauge coupling, rrig is the gluino mass, m~ h and m 4 ~ are the 
sbottom and stop masses, and A t is the top trilinear (scalar) coupling. 

The logarithmic corrections to yb are positive, which leaves the finite corrections 
to provide for the correct overall negative Syb in order to realize Yukawa unification. 
The gluino contribution (in Eq.(7)) is positive for \x > for the case of same sign 
gaugino SBB mass term. Thus, the chargino contribution (in Eq.(7)) must play 
an essential role to provide the required negative contribution to Syb- This can be 
achieved with suitably large values of both and A t , which is the reason behind 
the requirement m 16 > 6TeV and A /mi G ~ —2.6 in the 5*0(10) model discussed in 
[3,4]. 

One can improve the situation immensely by considering the case of opposite 
sign gaugino soft terms which is allowed by the 4-2-2 model. Reference [8] shows 
the parameter space, corresponding to /i < 0,M 2 < 0, M 3 > 0, that gives Yukawa 
coupling unification with a sub-TeV sparticle spectrum, and which is consistent with 
all known experimental bounds including A(g — 2) /1 . Another possibility is to consider 
/i > 0, M 2 > 0, M 3 < 0. This choice also improves the situation about t-b-r Yukawa 
unification, since it provides a negative gluino loop contribution to Syb- 
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4 t-b-r Yukawa unification with ra# = ra# 



After scanning the parameter space listed in Eq.(5), and employing the GUT scale 
boundary condition for the gauginos in Eq.(6), we find that REWSB can occur, 
compatible with t-b-r Yukawa unification, even if we set = m# d . Note that this 
is very difficult to achieve if gaugino mass universality is imposed at Mqut [28]. In 
order to explain our findings, consider the following one loop renormalization group 
equations for m 2 Hu , m 2 Hd and the soft trilinear terms At t h,r [29]: 

16vr 2 ^m^ = 3X t -6( ?2 2 |M 2 | 2 -^ 1 2 |M 1 | 2 + ^ ) (8) 
167r 2 ^m^ = ZX b + X T -6&\M 2 \ 2 -^g*\Mi\ 2 -^g*S. (9) 



9 d 
16vr 2 -A 
dt 



I67i 2 -A b 
dt 



lQn 2 ^A T 
dt 

Here 

X t 

x b 
x T 
s 



A t 

+y t 

A b 

+Vb 
A T 



16 



-9 3 



3gl 



32 



9 2 3 M 3 



13 

~ 15 
26 

6g 2 2 M 2 + ^g\M x 



9i 



2A b y b y t 



18y 2 b + y 2 + y 2 



16 



-9s 



3gl 



15 



2A t y t y b + 2A T y T y b 



32 14 

—g 2 3 M 3 + 6g 2 2 M 2 + 



(10) 



11} 



Wr + 3y 2 b 



3g 2 2 
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+ 6A b y b y T + y T 



18 

6<7 2 2 M 2 + -g\M x . (12) 



2\y t \ 2 (m 2 



Hv +m^ 3 +mio)+2\A t \ 



%IV! 



2\y T \~[m 
TrfFj-m! 



+ m% + m 2 3c ) + 2\A b 



,2 



Hd + m ^ +m ^) +2 \A T \ 



m 



m 



h + Tr[m 2 



m-. 



(13) 
(14) 
(15) 

2m? c + m 2 ~ c +m? c ], (16) 



and Q 3 , u 3 , d^, L 3 , e\ denote the third generation squarks and sleptons. Also, g^ and 
Mj (z = 1,2,3) denote the gauge couplings and gaugino masses for U(l)y, SU(2) L 
and SU(3)c, and yj, Aj (j = t,b,r) are the third family Yukawa couplings and 
trilinear scalar SSB couplings, respectively. 

In Eqs. (8) and (9), terms with negative coefficients drive m 2 Hu and m 2 ^ to larger 
values in the evolution from Mqut to Mz, while the terms with positive signs clearly 
have the opposite effect. On the other hand, at scale Mz ,the electroweak symmetry 
breaking minimization condition at tree level requires that 



m 



m 2 H tan 2 /3 



M 2 



tan 2 (3 



(17) 



For moderate to large values for tan/3 and assuming that \m 2 H \ ^> Mf, we can 
get /x 2 ~ —Tn 2 H . In order to have REWSB, m 2 H must be driven to negative values. 
In the literature, two alternative ways are used to achieve rrijj < 0. In one case 
with y t > i/b, y T , fn 2 H ls driven to negative values at Mz by y t , as seen from Eqs.(8) 
and (13). In this case it is possible to find solutions for which = at 

M GU t- In a model with universal gaugino masses at M GU t, we need a very large 
m 16 > 6 TeV and and even larger trilinear SSB coupling, A ~ — 2.6m 16 , in order 
that y t (M GVT ) = y b (M GVT ) = y T (M GVT ). As we know, one has y t > y b > y T all 
the way from Mqut to Mz, and the SSB terms in Eqs.(13) and (14) renormalize 
more or less equally. This means that X t > X b , but the additional contribution 
from X T makes the contributions from the positive terms in Eq.(9) larger than the 
corresponding positive contributions in Eq. (8). As a result, m 2 H is driven to negative 
values earlier than m 2 H at Mz, which is phenomenologically not acceptable. To avoid 
this problem, one can impose m# u < m# d at Mqut- As we mentioned in Section 
3, with opposite signs for M2 and M3 at Mqut, t-b-r Yukawa unification becomes 
easier by appropriately choosing the sign for /z in Eq. (7) [4, 8], and one can have 
t-b-r Yukawa unification even for m 16 ~ 1 — 2 TeV. Notice that the relative sign of 
M2 and M3 affects the RG running of the trilinear SSB couplings in Eqs.(10)-(12). 
In most cases, we find A t ~ A b , while A b and A T have opposite sign at Mz- We see 
from Figure la that A T is driven towards zero and then continues its evolution. This 
tendency helps in changing the direction of evolution for m# . 

As we will show in the next section, for precise t-b-r Yukawa coupling unification, 
rriH u = rriH d = m w < 1 TeV at M GUT . Hence, we can conclude that the SSB terms 
in Eqs.(13) and (14) are comparable. Since the values of the SSB terms for opposite 
signs of gaugino masses are around 1 TeV, the contribution from X T in Eq. (9) is not 
very significant. The coupling y t is slightly larger than y b during RG running, but in 
our case it is enough to drive m 2 H negative at the weak scale. 

The interplay among these various contributions is neatly summarized in Figure 
1. In the left panel we plot the RG evolution of the three soft trilinear couplings, 
while the right panel shows the corresponding evolution of vr? H and m 2 H . In order 
to have phenomenologically acceptable EWSB we need to have m 2 H more negative 
than rn 2 Hd at EWSB scale. The relation at high scale between m 2 Hu and m 2 H is not 
very important. Figure lb shows that at the EWSB scale m 2 H is more negative than 
m 2 H . In the text we use Mz as the scale of EWSB. 
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Figure 1: Evaluation of trilinear SSB couplings (left panel) and the MSSM Higgs 
(right panel) SSB mass terms corresponding to the benchmark point 1 in Table 1. 

5 Higgs Boson Mass and Sparticle Spectroscopy 

In order to quantify Yukawa coupling unification, we define the quantity R as, 

max{y t ,y b ,y T ) 

K = — r (18) 

nun(2fe, y b , y T ) 

In Figure 2 we show the results in the R — mi 6 , R — Myz, R — tan/3, R — A /miQ, 
R — mio and M1/2 — mie planes. The panels correspond to the following choice of 
parameters: fi > 0, M 2 > and M3 : M 2 : Mi = —2 : 3 : 1. The gray points are 
consistent with REWSB and neutralino LSP. The green points satisfy the particle 
mass bounds and constraints from BR(B S — > BR{b — > 57) and BR(B U — > 

rv T ). In addition, the green points do no worse than the SM in terms of (g — 2)^. The 
brown points belong to a subset of the green points and satisfy the WMAP bounds 
on neutralino dark matter abundance. 

In the R - m\§ plane of Figure 2 we see that we can realize perfect Yukawa 
unification consistent with all constraints mentioned in Section 2. This is possible 
because we can implement Yukawa unification for relatively small values of m 16 (~ 
1 TeV). This is more than an order of magnitude reduction in the mi6 values required 
for Yukawa unification, compared with the case /1 > and universal gaugino masses. 
In the R — Mi/ 2 plane of Figure 2, we see that employing the boundary conditions 
for gauginos presented in Eq. (6), perfect t-b-r Yukawa unification prefers heavier 
(> 1.2 TeV) values for Mi/ 2 . We can also predict that tan/3 m 47. In this case for 
Yukawa unification up to 5% level, 45 < tan/3 < 48. 

In the R — A /miQ plane of Figure 2 we see that in the present 5*0(10) case 
with 10% or better t-b-r Yukawa unification, we obtain a relaxation of Aq values 
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mio(TeV) mi 6 (TeV) 



Figure 2: Plots in R — rriiQ, R — M1/2, R — tan/3, R — A /m w , R — mi and Mi/ 2 — mi 6 
planes. Panels correspond to the following choice of parameters: fi > 0, M 2 > 
and M 3 : M 2 : M x = -2 : 3 : 1. Gray points are consistent with REWSB 
and neutralino LSP. Green points satisfy particle mass bounds and constraints from 
BR(B S — > BR{b — > sj) and BR{B U — > tu t ). In addition, we require that 

green points do no worse than the SM in terms of (g — 2) M . Brown points belong to 
a subset of green points and satisfy the WMAP bounds on neutralino dark matter 
abundance. 
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Figure 3: Plot in R — nih plane. The color coding is the same as that used in Fig. 2. 
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Figure 4: Plots in niq — nig, mf — m^o planes. The color coding is the same as that 
used in Fig. 2. In addition, yellow points represent Yukawa unification within to 
10%. Brown points in trig — m g plane that are compatible with WMAP bound on 
relic abundance are a subset of yellow points in nig — nig plane. 



similar to the SU(A) C x SU(2) L x SU(2) R model in [8], with -2.5 < A /m 16 < 1 for 
/i < 0, and —1.5 < A /m ie < 3 for fi > 0. The R — m w plane shows that precise 
unification prefers 0.2 TeV < m±o < 1 TeV. We see also that relaxing the Yukawa 
unification condition relaxes the upper bound for mio- In order to better visualize 
the magnitude of the sparticle masses consistent with t-b-r Yukawa unification, we 
present our results in the Mi/ 2 — m 16 plane, where the yellow points correspond to 
Yukawa unification better than 10%. 
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In Figure 3 we show the results in the R — rrih plane. It is most interesting to 
note that demanding precise t-b-r Yukawa unification for ix > 0, M 2 > and M 3 < 
allows us to predict the light CP even Higgs mass in a very narrow interval, 122 GeV 
< rrih < 124 GeV. For Yukawa unification of up to 5%, it becomes 120 GeV < rrih < 
125 GeV. We expect a theoretical uncertainty of about 3 GeV in the calculation of 
the light Higgs mass. 

In Figure 4 the rrig — rrig panel shows that t-b-r Yukawa unification in our scenario 
predicts masses for the gluino and the first two family squarks which lie somewhat 
beyond the current ATLAS [30] and CMS [31] bounds. 

In the present framework, the WMAP constraint on the relic dark matter abun- 
dance is only satisfied by the neutralino-stau co-annihilation scenario. We know that 
in neutralino-stau coannihilation, the neutralino and stau masses are degenerate to 
a good approximation. Instead of introducing a new color in the figure, we decided 
just to draw the unit line, and then it is easy to understand that points close to this 
line will give good relic abundance through neutralino-stau coannihilation. From the 
plot in m? — myX) plane of Figure 4, we see that the LSP neutralino mass is greater 
than or of order 300 GeV. Other co-annihilation scenarios may emerge after doing a 
more thorough analysis, but this remains to be seen. 

Finally, in Table 1 we present some benchmark points with fi > 0, M 2 > and 
M3 < 0. All of these benchmark points satisfy the various constraints, except possibly 
for a small discrepancy with the WMAP bound on relic dark density. Points 1 and 
2 depict solutions corresponding to minimum R and best VL CDM h 2 values. Point 2 
describes a solution with m 16 = m w . Point 3 displays a solution with a Higgs boson 
mass close to 124 GeV. The theoretical uncertainty in this estimate is 2-3 GeV. 

6 Conclusion 

We have reconsidered t-b-r Yukawa unification within a slightly revised framework 
in this paper. The main difference from most previous investigations stems from the 
assumptions we make related to the soft supersymmetry breaking parameters. First, 
we assume that the MSSM gauginos have non-universal masses, which are related to 
one another at M G ut by some appropriate SO(10) group theory factors. Second, we 
set the masses of the two MSSM Higgs doublets to be equal at Mqut- Overall, this 
means that we effectively have one parameter less than in the standard approach to 
SO (10) Yukawa unification. 

The ramifications of these slightly different assumptions for TeV scale physics 
turn out to be quite startling, with the low energy predictions being very different 
from previous studies. We find, for instance, that t-b-r Yukawa unification solutions 
at 5% level or better exist in our model with > and M 2 > 0. 

These solutions, obtained using the ISAJET software, are compatible with all 
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122.3 


122.1 


124 


mu 


408 7 




OOc/.O 


m A 


4U0 




000. o 


vn tta- 


420.5 


404.7 


568.3 


X-1,2 

TO-^o 

A3, 4 

m xf, 2 
m~ g 


609 3 2227 3 


^46 7 1 930 4 


8S3 4 3487 6 


2230.4, 3314.7 


1933.7, 2986.2 


3490.7, 4503.2 


ZZ0U.4, oZol.l 


1 Q^Q 8 OQ^fi 


oOoo.O, 4405.0 


r Jin (I 

5413.9 


4919.5 


7217.9 




5251.4, 4766.1 


4895.8, 4375.4 


7176.2, 6386 




3839. 2, 4665.7 


3453.1, 4211.6 


4967.6, 6120.1 


TO J 

d L,R 

mr 

01,2 


5352, 4759.7 


4896.5, 4369.9 


7176.7, 6378.7 


3969.3, 4640.1 


3589.8, 4188 


5168.1, 6088.1 


TO^ 


2784.8 


2585 


3844.3 


TOi> 3 


2639.9 


2421.7 


3605.8 


m e L ,R 


2789.5, 1355.5 


2589.3, 1359.5 


3849.1, 2016.8 


m fl 2 


609.4, 2650.7 


548.9, 2431.6 


842.4, 3621.2 


A(<? - 2)„ 


0.68 x 10- 10 


0.84 x 10- 10 


0.35 x 10- 10 


0"s/(pb) 


0.37 x 10~ 9 


0.59 x 10~ 9 


0.46 x 10~ 10 


o"5D (pb) 

^CDMh 2 


0.47 x 10~ 8 


0.89 x 10~ 8 


0.63 x 10~ 9 


0.176 


0.189 


0.55 


R 


1.01 


1.02 


1.01 



Table 1: Sparticle and Higgs masses (in GeV) for fi > 0. Points 1 and 2 depict 
solutions corresponding to minimum R and best VLcDNih 2 values. Point 2 corresponds 
to the case when m 16 = m 10 . Point 3 displays a solution with the heaviest mass for 
CP even light Higgs boson compatible with Yukawa unification better than 1%. 
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experimental observations, as well as the WMAP dark matter constraint (through 
stau co- annihilation). The masses of the gluino and first two family squarks are found 
to lie in the 2.7-5 TeV range, while the lightest stop (top squark) weighs at least 2 
TeV or so. Finally, with 5% or better Yukawa unification, the lightest Higgs boson is 
predicted to have a mass of around 120 - 125 GeV mass range, (with an uncertainty 
of ±3 GeV). The MSSM parameter tan/3 is around 45 - 47. 
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